Infrared optical response of strongly correlated cuprates: 
the effects of topological phase separation. 
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Cuprates are believed to be unconventional systems which are unstable with regard to a self- 
trapping of the low-energy charge transfer excitons with a nucleation of electron-hole (EH) droplets 
being actually the system of coupled electron CuO^ - and hole CuO^ - centers having been glued 
in lattice due to a strong electron- lattice polarization effects. The hole/electron doping into parent 
cuprate is likely to be a driving force for a growth of primary EH droplets with a gradual stabilization 
of a single, or multi-center skyrmion-like EH Bose liquid collective mode, followed by a first order 
phase transition from parent insulating state into unconventional topological EH Bose liquid phase 
(EHBL). Nanoscopic electron inhomogeneity is believed to be inherent property of doped cuprates 
throughout the phase diagram beginning from EH droplets in insulating parent system and ending 
by a topological phase separation in EHBL phase. 

We examine the effects of electron inhomogeneity accompanying such a phase separation on IR 
optical conductivity. A simple model of metal-insulator composite and effective medium theory has 
been used to describe the static phase separation effects. The low-frequency dynamics of topological 
EHBL phase in a random potential in underdoped regime has been discussed in a quasiparticle 
approximation in frames of the memory function formalism. 

The effects of static and dynamic nanoscopic phase separation are believed to describe the main 
peculiarities of the optical response of doped cuprates in a wide spectral range. 

Keywords: cuprates, phase separation, optical conductivity, memory function 
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I. INTRODUCTION 

The nature of the optical response of strongly corre- 
lated cuprates is presently still a matter of great contro- 
versy. The unconventional behavior of cuprates strongly 
differs from that of ordinary metals and merely resembles 
that of a doped semiconductor. Moreover, the history of 
high T c 's itself shows that we deal with a transforma- 
tion of particularly insulating state in which the electron 
correlations govern the essential physics. New principles 
must be developed to treat these systems with its non- 
Fermi-liquid behavior. First of all we have to change the 
current paradigm of the metal-to-insulator (MI) transi- 
tion to that of an insulator-to-metal (IM) phase transi- 
tion. These two approaches imply essentially different 
starting points: the former starts from a rather simple 
metal scenario with inclusion of correlation effects, while 
the latter does from strongly correlated atomic-like sce- 
nario with the inclusion of a charge transfer. 

The copper oxides start out life as insulators in con- 
trast with BCS superconductors being conventional met- 
als. It is impossible to understand the behavior of the 
doped cuprates and, in particular, the origin of HTSC 
unless the nature of the doped insulating state is incor- 
porated into the theory. In particular, the Fermi liquid 
theory of the normal state and the BCS theory of the 
superconducting state, which are so successful for con- 
ventional metals, were not designed for doped insulators, 
and they do not apply to the HTSC. Consequently it is 
necessary to develop a new mechanism and many-body 
theory of HTSC. 

The problem of a doped insulator is sure much more 
complicated than it is implied in oversimplified ap- 



proaches such as an effective t- J- model. First, we deal 
with effects of electron inhomogeneity and phase separa- 
tion. These are currently considered to be inherent prop- 
erty of lightly doped or so-called under-doped cuprates. 
However, recent scanning electron micrography measure- 
ments clearly evidence the inherent electron inhomogene- 
ity in optimally doped BiaSraCaCu-jOgia" 

The emergence of electron inhomogeneity leads to a 
strong complication of the optical response and its anal- 
ysis due to a number of specific effects typical for opti- 
cally inhomogeneous systems. The discussion and mod- 
elling of such effects would be the main goal of the paper. 
Its structure is organized as follows. In Sec. II we dis- 
cuss the origin and evolution of electron inhomogeneity 
in doped cuprates. In Sec. Ill the effective medium theory 
is applied to describe the optical response of " quenched" 
electronic inhomogeneity. The effects of the dynamics of 
multi-center topological defect in a random potential are 
addressed in Sec. IV. 



II. CONDENSATION OF CT EXCITONS AND 
PHASE SEPARATION IN DOPED CUPRATES 

In case of cuprates we deal with systems which con- 
ventional ground state seems to be unstable with regard 
the transformation into a new phase state with a variety 
of unusual properties. It is worth noting the text-book 
example of BaBiC>3 system where we unexpectedly deal 
with the disproportionated Ba 3+ + Ba 5+ ground state in- 
stead of the conventional lattice of Ba 4+ cations. The 
bismuthate situation can be viewed as a result of conden- 
sation of charge transfer (CT) excitons, in other words, 
the spontaneous generation of self-trapped CT excitons 
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in the ground state with a proper transformation of lat- 
tice parameters. These problems are closely related with 
the hidden multistability intrinsic to each solid£ If the 
ground state of a solid is pseudo-degenerate, being com- 
posed of true and false ground states with each structural 
and electronic orders different from others, one might call 
it multi-stable. 

In our view, the physics of the doped cuprates, in- 
cluding superconductivity, is driven by a self-trapping 
of the CT excitons?^ Such excitons are the result of 
self-consistent charge transfer and lattice distortion with 
appearance of the "negative-?/" effect. In contrast with 
BaBi03 system where we deal with a spontaneous gen- 
eration of self-trapped CT excitons (STE) in the ground 
state, the parent insulating cuprates are believed to be 
near excitonic instability when the self-trapped CT exci- 
tons form the candidate relaxed excited states to struggle 
with the conventional ground stated In other words, the 
lattice relaxed CT excited state should be treated on an 
equal footing with the ground state. If the interaction 
between STE were attractive and so large that the cohe- 
sive energy W\ per one STE exceeds the energy En of 
one STE, the STE's and/or its clusters will be sponta- 
neously generated everywhere without any optical exci- 
tation, and be condensed to form a new electronic state 
on a new lattice structured 



A. Electron-hole droplets in cuprates: 
pseudo-impurity and phase separation regimes 

Cuprates are believed to be unconventional systems 
which are unstable with regard to a self-trapping of the 
low-energy charge transfer excitonsS^S with a nucleation 
of EH droplets being actually the system of coupled elec- 
tron CuO^" and hole CuO|~ centers having been glued 
in lattice due to a strong electron-lattice polarization ef- 
fects. Phase transition to novel hypothetically metal- 
lic state could be realized due to a mechanism famil- 
iar to semiconductors with filled bands such as Ge and 
Si where given certain conditions one observes a for- 
mation of metallic EH-liquid as a result of the exciton 
decay. 9 However, the system of strongly correlated elec- 
tron CuC>4~ and hole CuO^ - centers appears to be equiv- 
alent to an electron-hole Bose-liquid (EHBL) in contrast 
with the electron-hole Fermi-liquid in conventional semi- 
conductors. A simple model description of such a liquid 
implies a system of local singlet bosons with a charge of 
q = 2e moving in a lattice formed by hole centers^ Lo- 
cal boson in our scenario represents the electron counter- 
part of Zhang-Rice singlet, or two-electron configuration 
b\ g 1 Aig. Naturally, that conventional electron CuOj" 
center represents a relaxed state of a composite system: 
"hole CuO^ - center plus local singlet boson", while the 
"non-retarded" scenario of a novel phase is assumed to 
incorporate the unconventional states of electron CuOj - 
center up to its orbital degeneracy. 



Homogeneous nucleation implies the spontaneous for- 
mation of EH droplets due to the thermodynamic fluc- 
tuations in exciton gas. Generally speaking, such a state 
with a nonzero volume fraction of EH droplets and the 
spontaneous breaking of translational symmetry can be 
stable in nominally pure insulating crystal. However, 
the level of intrinsic non-stoihiometry in oxides is sig- 
nificant (one charged defect every 100-1000 molecular 
units is common). The charged defect produces ran- 
dom electric field, which can be very large (up to 10 8 
Vcm -1 ) thus promoting the condensation of CT exci- 
tons and the nucleation of EH droplets. Deviation from 
the neutrality of the Cu02 layers implies the existence 
of additional electron, or hole centers that can be the 
natural centers for the inhomogeneous nucleation of the 
EH droplets. Such droplets are believed to provide the 
more effective screening of the electrostatic repulsion for 
additional electron/hole centers, than the parent insu- 
lating phase. As a result, the electron/hole injection to 
the insulating cuprate due to a nonisovalent substitu- 
tion as in La2- a; Sr 2 ;Cu04, Nd2- a; Ce a ;Cu04, or change in 
oxygen stoihiometry as in YBa2Cu30g+a;7 La2Cu04_a, 
La2Cui_ a: Li 2 ;04, or field-effect is believed to shift the 
phase equilibrium from the insulating state to the un- 
conventional electron-hole Bose liquid, or in other words 
induce the insulator-to-EHBL phase transition. 

The optimal way to the nucleation of EH droplets in 
parent system like La 2 Cu04, YBa2Cu306 is to create 
charge inhomogeneity by nonisovalent chemical substitu- 
tion in Cu02 planes or in "out-of-plane stuff" , includ- 
ing interstitial atoms and vacancies. This process re- 
sults in an increase of the energy of the parent phase and 
creates proper conditions for its competing with others 
phases capable to provide an effective screening of the 
charge inhomogeneity potential. The strongly degener- 
ate system of electron and hole centers in EH droplet 
is one of the most preferable ones for this purpose. At 
the beginning (nucleation regime) an EH droplet nucle- 
ates as a nanoscopic cluster composed of several number 
of neighboring electron and hole centers pinned in the 
CuC>2 plane by disorder potential. Of course, the condi- 
tions for such a nucleation are delicately determined by 
the crystal and chemical surroundings of the Cu02 planes 
("out-of-plane stuff"). Both the electron and hole cen- 
ters at the atomic-molecular level and the EH droplets 
at the nanoscopic scale level are the elementary building 
blocks for an unconventional electron-crystalline struc- 
ture of the cuprates with the nonisovalent substitution. 
EH droplet is a nanoscopic center for a pinning of the var- 
ious quasi-local hybrid charge-spin-vibronic modes with 
a definite relaxation time and spatially non-uniform DOS 
distribution or even internal micro-phase separation. 

An occurrence of the EH droplets as the domains of 
the novel phase manifests itself in various properties 
of the cuprates with a non-isovalent substitution even 
at the low content of the nucleation centers. On the 
one hand, main features of this pseudo-impurity regime 
are determined by the intrinsic properties of the EH 



3 



droplets with expected derealization and metallic-like 
behavior at T > Tcoi localization with charge ordering 
at T < Tcoi an d other phase transformations up to a 
quasi-local bose-condensation at T < Tbs with an ap- 
pearance of the superconducting order parameter^! On 
the other hand, the real properties of the system will be 
determined by the peculiar geometrical factors as a rel- 
ative volume fraction of the EH droplets, their average 
size and appropriate dispersion, and, possibly, by their 
form and/or texture. Numerous examples of the uncon- 
ventional behavior of the cuprates in the pseudo- impurity 
regime could be easily explained with taking into account 
the inter-phase boundary effects (coercitivity, the mo- 
bility threshold, oscillations, relaxation etc.) and corre- 
sponding characteristic quantities. In this connection one 
should emphasize an occurrence of the characteristic tem- 
peratures Tb for the start of the boundaries motion, or for 
the loss of the stability of either geometrical droplet con- 
figuration. In general, such phenomena should accom- 
pany any partial phase transition for a phase-separated 
system, that could result in an appearance of a peculiar 
doublet structure of the temperature anomalies. 

EH droplets can manifest itself remarkably in various 
properties of the cuprates even at small volume fraction, 
or in a "pseudo-impurity regime" . Main features of this 
"pseudo-impurity regime" are determined both by the 
intrinsic properties of the EH droplets and by the pecu- 
liar geometrical factors as a relative volume fraction of 
droplets, their average size and appropriate dispersion, 
and, possibly, by their form and/or texture. Numerous 
examples of the unconventional behavior of the cuprates 
in the pseudo-impurity regime could be easily explained 
with taking into account the inter-phase boundary effects 
(coercitivity, the mobility threshold, oscillations, relax- 
ation etc.) and corresponding characteristic quantities. 

Under increasing doping the "pseudo-impurity regime" 
with a relatively small volume fraction of EH droplets can 
gradually transform into a micro- (electronic) or macro- 
chemical) "phase-separation regime" with a sizeable vol- 
ume fraction of EH droplets. 

The above viewpoint on the phase separation (PS) phe- 
nomena in cuprates is, in general, compatible with the 
pioneer ideas by Emery and Kivelsoni2iiii4 and some 
other model approaches. So on the basis of the neutron 
scattering data Egamii^ conjectured an appearance of a 
nano-scale heterogeneous structure which is composed of 
a plenty mobile-carrier existing region of metallic con- 
ductivity and semi-localized scarce carrier region with 
antiferromagnetic spin ordering. 

Phase separation is now widely discussed as an impor- 
tant phenomenon accompanying the high-T c supercon- 
ductivity. The pseudo-impurity regime can be termed as 
a nanoscopic phase separation. Macroscopic, or chemical 
phase separation with the x-rays and neutrons detectable 
domains of a new phase exactly takes place in super- 
oxygenated La2Cu04+^i&, due to the enhanced mobil- 
ity of negatively charged oxygens which can compensate 
charge unbalance introduced by the hole separation. A 



tendency towards such a phase separation has also been 
repeatedly observed in other cuprates. 

At once, the well developed inhomogeneity and phase 
separation are inconsistent with the excellent coherence 
of the lattice demonstrated by the sharpness of many 
of the Bragg peaks (for example, coherence length over 
1000A in La 2 _ 2; M :c Cu0 4 ii). However, an appearance 
of the clearly seen nonphonon X-ray diffuse scattering 
with unconventional temperature and q-dependence 1 ? 
evidences an occurrence of the local or quasi-local 
(nanoscale) static and/or dynamic charge inhomogene- 
ity as a possible result of the electronic phase separation. 
Lack of significant elastic neutron diffuse scatteringii is 
compatible with the mainly electronic nature of the in- 
homogeneity in doped cuprates. 



B. Topological phase separation in doped cuprates 

The phase separation and electron inhomogeneity are 
currently considered to be inherent property of lightly 
doped or so-called under-doped cuprates. Our scenario 
implies the existence of a critical doping level for the 
first-order phase transition from parent insulating phase 
to EHBL phase. EHBL represents the phase state with 
a condensation of CT excitons. Without any doping 
such a phase is equivalent to a half-filled hard-core Bose- 
Hubbard (hc-BH) system. Does new phase is homoge- 
neous or not? We argue that EHBL phase is not only 
inherently inhomogeneous but this inhomogeneity gov- 
erns the low-energy physics of doped cuprates. Indeed, 
the model of the EH Bose-liquid implies that the doped 
cuprates are proved to be in the universality class of the 
pseudo-spin 2D systems which description incorporates 
the inhomogeneity centers like topological defects to be 
a natural element of essential physics. In addition, it 
should be noted that the self-trapping of CT excitons 
with nonzero electric dipole moment is characterized by 
the nucleation of vortex-like configurations providing the 
minimum of the dipole-dipole coupling. In other words, 
the nucleation of the EH Bose liquid implies the forma- 
tion of different topological textures. However, below we 
concentrate on the deviation from the half-filling to be 
seemingly the main driving force of the topological inho- 
mogeneity in doped cuprates. 

The hc-BH- model (quantum lattice Bose gas) has a 
long history and has been suggested initially for conven- 
tional superconductors^, quantum crystals such as A He 
where superfluidity coexists with crystalline order^Si 
Afterwards, the Bose-Hubbard model has been studied 
as a model of the superconductor-insulator transition in 
materials with local bosons, bipolarons, or preformed 
Cooper pairs iSLSS The most recent interest to the system 
of hard-core bosons comes from the delightful results on 
Bose-Einstein (BE) condensed atomic systems. 

One of the fundamental hot debated problems in 
bosonic physics relates the evolution of the charge or- 
dered (CO) ground state of 2D hard-core BH model with 
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doping away from half-filling. Numerous model studies 
steadily confirmed the emergence of " supersolid" phases 
with simultaneous diagonal and off-diagonal long range 
order while Penrose and Onsager— were the first showing 
as early as 1956 that supersolid phases do not occur. 

The most recent quantum Monte-Carlo 
simulations2^25i2& found two significant features of 
the 2D Bose-Hubbard model with a screened Coulomb 
repulsion: the absence of supersolid phase at half-filling, 
and a growing tendency to phase separation (CO-BS) 
upon doping away from half-filling. Moreover, the 
checkerboard supersolid phase appears to be unstable 
thermodynamicallyi^i2^2& Batrouni and Scalettar stud- 
ied quantum phase transitions in the ground state of 
the 2D hard-core boson Hubbard Hamiltonian and have 
shown numerically that, contrary to the generally held 
belief, the most commonly discussed "checkerboard" 
supersolid is thcrmodynamically unstable and phase 
separates into solid and superfluid phases. The physics 
of the CO-BS phase separation in Bose-Hubbard model 
is associated with a rapid increase of the energy of a 
homogeneous CO state with doping away from half- 
filling due to a large "pseudo-spin-flip" energy cost. 
Hence, it appears to be energetically more favorable to 
"extract" extra bosons (holes) from the CO state and 
arrange them into finite clusters with a relatively small 
number of particles. Such a droplet scenario is believed 
to minimize the long-range Coulomb repulsion. 

What is, however, the detailed structure of the CO+BS 
phase separated state? In the paper— a topological sce- 
nario of CO+BS phase separation has been proposed. 
The extra bosons/holes doped to a checkerboard CO 
phase of 2D boson system are believed to be confined 
in the ring-shaped domain wall of the skyrmion-likc 
topological defect which looks like a bubble domain in 
antifcrromagnet. 2 - 8 This allows us to propose the mecha- 
nism of 2D topological ( CO+BS) phase separation when 
the doping of the bare checkerboard CO phase results in 
the formation of a multi-center topological defect, which 
simplified pseudo-spin pattern looks like a system of bub- 
ble CO domains with Bose superfluid confined in charged 
ring-shaped domain walls. 

Such a skyrmion scenario in hc-BH model allows us to 
draw several important conclusions. The parent checker- 
board CO phase gradually looses its stability under bo- 
son/hole doping, while a novel topological self-organized 
texture is believed to become stable. In other words, 
away from half-filling one may anticipate the nucleation 
of topological defect in the unconventional form of the 
multi-center skyrmion-like object with ring-shaped Bose 
superfluid regions positioned in antiphase domain wall 
separating the CO core and CO outside of the single 
skyrmion. Specific spatial separation of BS and CO or- 
der parameters that avoid each other reflects the compe- 
tition of kinetic and potential energy. Such a topological 
( CO+BS) phase separation is believed to provide a mini- 
mization of the total energy as compared with its uniform 
supersolid counterpart. 



The most probable possibility is that every domain 
wall accumulates single boson, or boson hole. Then, 
the number of centers in a multi-center skyrmion nucle- 
ated with doping has to be equal to the number of doped 
bosons/holes. 

Apparently, the boson/hole doping is likely to be a 
driving force for a nucleation of a single, or multi-center 
skyrmion-like self-organized collective mode, which may 
be (not strictly correctly) referred to as multi-skyrmion 
system akin in a quantum Hall ferromagnetic state of 
a two-dimensional electron gas 2 ^. We may characterize 
an individual skyrmion by its position (i.e., the center 
of skyrmionic texture), its size (i.e., the radius of do- 
main wall), and the orientation of the in-plane compo- 
nents of pseudo-spin (U(l) degree of freedom). An iso- 
lated skyrmion is described by the inhomogeneous dis- 
tribution of the CO parameter, or staggered boson den- 
sity, charge order parameter characterizing the deviation 
from the half-filling, and superfluid order parameter with 
a modulus and phase. As regards the charge distribution, 
the multi-center topological defects has much in common 
with a charge-density wave, Wigner crystal or Skyrmion 
crystal. 

It is worth noting that multi-center skyrmions, or 
multi-skyrmion topological defects in non-linear er-model 
were regarded by Belavin and Polyakov. 7 Interestingly, 
that such entities may be considered as a system of non- 
interacting quasiparticles. 

It seems likely that any doped particle (boson/holes) 
results in a nucleation of a new skyrmion, hence its den- 
sity changes gradually with particle doping. Therefore, 
as long as the separation between skyrmionic centers is 
sufficiently large so that the inter-skyrmion interaction is 
negligible, the energy of the system per particle remains 
almost constant. This means that the chemical poten- 
tial of boson or hole remains unchanged with doping and 
hence apparently remains fixed. 

Depending on the doping level we may distinguish 
three types of topologically phase-separated EHBL in 
doped cuprates: i)"underdoped" EHBL with weakly cou- 
pled skyrmions, when their behavior is ruled mainly by 
random impurity potential; ii) "optimally doped" EHBL 
with strongly coupled skyrmions, forming skyrmion liq- 
uid, or crystal; iii) "overdoped" EHBL with a complex 
topological texture that hardly reduces if any to a system 
of well separated quasiparticles. The illustrative picture 
of different phase separation regimes in doped cuprates 
is presented in Fig^ 



III. IR OPTICAL RESPONSE OF THE DOPED 
CUPRATE. STATIC TOPOLOGICAL PHASE 
SEPARATION AND EFFECTIVE MEDIUM 
THEORY 

The hole/electron doping into parent cuprate is likely 
to be a driving force for a growth of primary electron-hole 
droplets with a stabilization of a single, or multi-center 
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FIG. 1: Illustrative picture of different phase separation regimes in doped cuprates 



skyrmion-like EHBL collective mode, followed by a first 
order phase transition from parent insulating state into 
unconventional topological EHBL phase. The appropri- 
ate evolution of electron structure should manifest itself 
in optical properties. We believe, that in a broad spec- 
tral range the basic features of the optical response of SC 
oxides can be consistently accounted for proceeding from 
the unified physical picture of their intrinsic electronic 
inhomogeneity. Below we focus on optical response in 
infrared. First we address the response of doped cuprate 
with a static, or "quenched" phase separation when the 
crude reasonable model is that of a granular composite, 
whose constituents have the complex dielectric functions 
Ei(u>), where i labels the species of the grains. What is to 
be stressed is that the resulting dielectric function e e ff 
of the composite does not reduce to the simple weighted 
average of the form pe\ + (1 — p) £2, where p is the vol- 
ume fraction of its first component. The correct averag- 
ing can be performed within the effective medium theory 
and takes into account the specific effects of the resonant 
scattering and absorption of light at the interfaces of the 
grains. The corresponding surface plasmon resonances 
in the optical spectra of the composite media are also 
referred to as the geometric ones, because their energy is 
governed by the shape of the inhomogeneities. 

Optical measurements in midinfrared (MIR) region 
provide an important information as regards the low- 
energy excitations in cuprates. The nature of MIR band 
in insulating cuprates remains to be one of the old mys- 
teries of the cuprate physics. 

At first glance, overall MIR band in insulating cuprates 
may be assigned to different dipole-allowed CT transi- 
tions in EH droplets which certain volume fraction is 
surely contained in the nominally pure parent cuprates. 
Generally speaking, the low-frequency optical response 
of EH droplet is composed from two main contributions, 
that of intra-band two-center local boson transitions and 
one-center 1 A\ — 1 E u CT transitions. The former can 



form a metallic-like band with the maximal width of the 
order of 0.1 eV, while the latter represents the complex 
optical response of single PJT center which typical en- 
ergy falls in the interval of several tenth of eV. However, 
the EH droplets are scattered in a highly polarizable sur- 
rounding that turns the calculation of their optical re- 
sponse into a very complicated problem. Typical size of 
the EH droplets seems to fall in the nanoscopic scale. 
Hence, when considering the phenomena with a relevant 
characteristic length, e.g. infrared optical absorption (A 
= 0.5 - 1.5 fim), doped cuprates can be roughly regarded 
as a binary granular medium, composed of two compo- 
nents. 

The idea of real-space inhomogeneity is most naturally 
expressed in terms of an effective medium theory (EMT) , 
which aims to evaluate some property of a composite, 
those of pure components being given. First, we would 
like to shortly overview the effective medium theory^ 
which appears to be a powerful tool for the quantitative 
description of the optical response of inhomogeneous sys- 
tems. In its simplest form, the EMT equation for effective 
dielectric function e e f / of the two-component inhomoge- 
neous system reads as follows^: 

/ dVp(e) V £ = £ ; S! r = 0, (1) 

J H £ eff +L t (e- S eff ) 

v '- 1 

where the integration runs over the volume of the sample, 
p(e) = pi8{e{r) — E\) + p2S(e(r) —£2). We consider 
the two components of binary composite on equal foot- 
ing with pi 2 being the volume fractions, £1,2 and Li the 
dielectric functions and the depolarization factors of the 
grains of its constituents, respectively. 

It is worth to emphasize that in effective medium ap- 
proximation the volume of the droplet does not enter 
the calculations, but is accounted implicitly in the valid- 
ity range of the theory, restricting the mean size of the 
droplet to be smaller than the wavelength. Hence, within 
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the EMT, the volume fractions p\^ can only change 
through the number of the droplets rather than due to 
the variation of their sizes. In general, pi.2 are deter- 
mined by thermodynamical conditions and depend on 
temperature, pressure, and other external factors. The 
approach^, we employed here to lend more plausibility to 
this physically transparent EMT scheme, is to take into 
account the natural difference between the core and the 
surface properties of the EH droplets using the standard 
expression for the polarizability of the coated ellipsoid. 32 
The spectra of nanoscopically disordered metal- 
insulator media can display some specific features due to 
so-called geometric resonances (surface plasmon, or Mie 
resonances) 3 ^, that have no counterpart in homogeneous 
systems. These arise as a result of resonant behaviour of 
local field corrections to the polarizability of the granu- 
lar composite and are governed to a considerable extent 
by the shape of the grains. The frequency of geomet- 
ric resonance is then easily obtained as the one at which 
the polarizability of small particle diverges. For the case 
of spherical metallic droplets embedded in the insulat- 
ing matrix with a constant dielectric permittivity Ed this 
leads to the equation: 



e(oj) m + 2e d = 0, 
whence the resonance frequency is 

w r = =. 
VI + 2e d ' 

if simple Drude's expression 

,2 



S = S() 



lu (lu + i7) 



(2) 



(3) 



(4) 



with the plasma frequency u p and the damping param- 
eter 7 is assumed for the metallic permittivity. In the 
case of arbitrary ellipsoid there are three different prin- 
cipal values of polarizability and generalizes to 



the in-plane semi-axes ratio denoted as a and that of 
out-of-plane and major in-plane semi-axes as (3. These 
geometrical parameters enter the theory through the de- 
polarization factors and were found to affect strongly the 
dielectric function of doped cuprate. 

Left hand side panel in Fig|21 clearly shows the evo- 
lution of the low-frequency optical response (MIR-band) 
given small volume fraction of metallic- like (u> p = 1.9 
eV) disc-shaped droplets with varying effective relaxation 
rate and volume fraction. For comparison we demon- 
strate the predictions of a simple additive two-fluid model 
a(uj) = pametai + (1 - p) o insulator (dotted curves). One 
may see that the effect of EH droplets dispersed in a po- 
larizable insulating matrix results in a crucial rearrange- 
ment of the optical response with the appearance of novel 
characteristic features such as Mie resonances. It is worth 
noting that in terms of EMT approach the appearance of 
metallic-like EH droplets in insulating matrix leads to the 
crucial red-shift of the effective insulating gap of doped 
cuprate. In a sense, for lightly doped cuprates one might 
say about two characteristic insulating gaps. 

The model in its simplest form— was found able to 
reproduce all essential features of the transmittance 35 , 
optical conductivity er(w), and EELS spectra— for 214 
system. It is easy concluded after simple comparison of 
the data presented in central and right hand side pan- 
els of Fig[3 Substantial difference in the spectral and 
doping dependence of optical conductivity for thin-film 35 
and bulk samples^ is easily explained assuming different 
shape of metallic and dielectric regions in both materials. 
New peaks in ct(cl>) and absorption spectra, that emerge 
in the mid-infrared range upon doping are attributed to 
surface plasmon resonances and its interference with dif- 
ferent charge transfer transitions in EH droplets. 

Overall, one may conclude that the MIR optical re- 
sponse of insulating cuprates strongly supports the EH 
droplet scenario with an intrinsic metallic-like transport. 



UJ,, 



Ed 



U(e d - 1)' 



= 1, 2, 3, (5) 



where Li are three shape-dependent depolarizaton fac- 
tors. The Drude approximation is reasonable and is more 
general that it may seem, since the form can be regarded 
as a limiting one of the universal expression for a dy- 
namical conductivity in terms of memory function with 
Up and 7 to be effective parameters. 

An effective medium theory has been applied for the 
simulation of the optical spectra of La2- x Sr :E Cu04 with 
x varying in a wide ranged (see Fig|2J). Imaginary part 
of dielectric function of parent insulating La2Cu04 was 
fitted to experiment— by the sum of three Gaussians in a 
broad spectral range to ensure the validity of its real part, 
derived via Kramers-Kronig transformation. The phonon 
bands have been neglected throughout the calculations. 
The droplets are assumed to have ellipsoidal shape with 



IV. IR OPTICAL RESPONSE OF THE DOPED 
CUPRATE. DYNAMICS OF TOPOLOGICAL 
DEFECT IN A RANDOM POTENTIAL 



The topological structure of EHBL must be considered 
as being largely dynamic in nature. Underdoped EHBL 
can be considered to form a system of weakly coupled 
quasiparticles (skyrmions) exposed to the action of a ran- 
dom potential. We believe, that the model of 2D Bosc 
gas in a random potential offers the most close physi- 
cal analogy to our case, and we adopt it as the effective 
model to describe the effects of dynamics of topological 
defect in the cuprates. 

Below we'll address the approximation of slow topo- 
logical dynamics with characteristic energies much lower 
than the effective insulating gap of doped cuprate. 



0.2 0.4 0.6 0.8 1 
Photon Energy (eV) 



12 3 4 

Photon Energy (eV) 



12 3 4 

Photon Energy (eV) 



FIG. 2: Optical conductivity of doped 214 cuprate in frames of the effective medium theory. Left hand side panel shows 
the evolution of the low-frequency optical response given small volume fraction of metallic-like droplets and varying effective 
relaxation rate. Central panel: calculated optical conductivity for different volume fractions of metallic phase. The numbers 
against the curves stand for the relative volume fraction p of the quasi-metal phase. Right hand side panel: measured optical 
conductivity of La2-;cSr x Cu04 for different x [34|]. The numbers against the curves stand for the doping level. 



A. The memory function technique. 

An adequate formalism for the description of the over- 
damped charge dynamics in presence of a strong dis- 
order was elaborated by W. Gotze and coworkers^ 6 * 3 ^, 
who showed, that the optical conductivity a(u>) at both 
sides of Anderson transition can be expressed in a general 
form: 



a' + i a" 



M{uj) 



(6) 



where lu p is the plasma frequency, M(oj) = M' + i M" is 
the memory function. With M = i 7 = const this ex- 
pression reduces to the simple Drude formula (??). Thus, 
the Drude model can be understood as a crudest approxi- 
mation, when the processes of the current dissipation act 
equally in the whole spectral range and do not depend 
on the frequency. 



The real part of the conductivity a', henceforth de- 
noted as a, may be parameterized in terms of the so 
called generalized Drude formula 36 , 



a(w) 



w 2 + r(w) 2 ' 



(7) 



where the renormalized plasma frequency £l p (u) and the 
relaxation rate T(uj) have the form: 



Q p (oj) = ^ p /\/^H, 7 (w) - M"(u)/ti(u), (8) 
r](w) = l + M'(u))/w. (9) 

The cornerstone result of the approach is based on the 
mode coupling self-consistency relation of the memory 
function and the impurity potential in the momentum 
spaced: 



M(lo) = 



$oO, u) + M(u>)) 



l + M(u)$o(k,u + M(w))/x(k = 0,uj) 
I 



k 2 (\U(q)\ 2 )dk, 



(10) 



where n, m are the concentration and the mass of the 
carriers, $0 is the density-density correlation function in 
absence of the impurities, expressed through the gener- 
alized susceptibility n(q,z): 



1 



®o(q,z) = - (n(q,z) - K(q) z =o). 



(11) 



Within the linear response theory n(q, z) is obtained from 



the familiar expression: 



lim > 



f{£k+ q ) - f(e q ) 

£k+q — £q — U — i 



(12) 



where is the energy of the particles and / is their 
partition function. Other quantities in Eq. <|10[) are the 
concentration n and the mass m of the charge carriers. 

Now, dwell upon the quantity (\U(q)\ 2 ). This is the 
average squared Fourier transform of the random poten- 
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tial U(r), defined a superposition of the partial impurity 
potential u, each centered at the sites iv 



U(r) = ^>(r- r *)' 



(13) 



where the sum is over the impurity sites. Hence, 

(\U(q)\ 2 ) = \u(q)\ 2 S (q) (14) 



»(«) = (£« 



£q (r ( - rj ) 



(15) 



where the brackets denote the averaging over the impu- 
rity sites. In what follows, we shall consider the case 
of random impurities. Then, the structure factor s(q) 
becomes trivial and reduces to rii, the concentration of 
impurities, whence 



(P(q)\ 2 ) = ni \u(q)f 



(16) 



In zero temperature limit, the machinery of the correla- 
tion functions of free 2D Bose gas is easily implemented 
because of the particularly simple analytic form of the 
generalized susceptibilitjs^: 



Ko(k, z) 



2ne k 



(17) 



where e k = k 2 /2m is the dispersion of free particles. 

To proceed further, it is convenient to introduce the 
dimensionless variables, labeled with tilda (~) to dis- 
tinguish them from their dimensional counterparts. We 
have chosen for the units of the mass and the charge 
the corresponding quantities for the free electron (m e , 
e). Aiming to simplify the calculations, we adopt the 
unit of the wave vector as q* — (Sim/ a*) 1 ^ , where 
a* = h 2 e/me 2 is the effective Bohr radius^, e is the static 
dielectric constant of the medium. The unit of energy is 
then naturally defined as E* = h 2 q* 2 /2m e . Furthermore, 
the complex optical conductivity takes the form: 



a + i a 



E* 



uj + M(u) ' 
1 n ZD e 2 
E* m 



(18) 
(19) 



where uj = lu/E*, M = M/E*, n 3 £> being the volume 
fraction of the droplets. Under moderate doping, we be- 
lieve this quantity to be directly related to the doping 
level x: 



n 3 D 



A 2 cf 



(20) 



where A is the lattice constant of the CuC>2 plane, d is 
the interplane spacing. 

Taking into account the numerical values of m e , e e 
and h, the final expressions for the adopted units read as 



follows: 

Q* 

E* 



3.62 • 10 8 x 



49.87 x 



22 ~\ V3 



me x 
eA 2 



e x 



eA 2 



2/3 



4.0 • 10 3 x 



/ e 2 e 2 x 
\ A 2 d 3 in 2 



eV, 



1/3 



(21) 
(22) 

n _1 cm _1 , (23) 



where we made use of the relation 1 eV 



1.824 



10 3 17 _1 cm~ 1 . Let us estimate their possible numerical 
values. Consider the situation, when the EHBL phase is 
composed of the heavy Bose particles carrying double ele- 
mentary charge, rh = 10, e = 2. Given the parameters of 
La2- 2; Sr a ,Cu04 unit cell, A ~ 4 A, d ~ 6 A, and the bare 
dielectric constant e = 40, we obtain for slight doping 
x = 0.01 the following units of energy and conductivity: 
E* ~ 294.5 cm- 1 , a* ~ 0.23- 10 3 ft" 1 cm" 1 . 

Making use of these notations, it is straightforward 
to write down the explicit expressions of the correlation 
functions of the theory: 



n(q,z) — [2n/E*) n(q, z) 
$(q,z) = (2n/E* 2 )$(q 7 ?) 
(z = E* z, q = q*q), 



(24) 
(25) 



where the frequency and wave vector dependencies are 
contained in dimensionless expressions: 



K(q,z) 



RPA 



§(q,z + M) 



q 4 + q-ti 2 ' 



q(z + M) 



(l + g 3 )(q + q 4 - (z + M)z) 



and 



n(k, z) 



RPA = 



n(k, z) 



1 + 2ne 2 K(k,z)/k 



(26) 



(27) 



(28) 



takes into account the 2D Coulomb interaction of quasi- 
particles in a random phase approximation. 

It is also convenient to isolate the dimensionless q- de- 
pendence in the impurity potential: 



«(«) = U f(q), 



(29) 



where Uq has the dimension of energy and / specifies 
the shape of the potential. Finally, the dimensionless 
analogue of Eq. I|10|) has the form: 



(30) 



Thus, the dimensionless memory function is governed by 
the only parameter, which is defined as: T = 4J7g n,/n, 
where Uq, measured in units of E* , is the strength of 
the potential (the depth of the potential well) , rij and n 
are the concentration of the impurities and the carriers, 
respectively. 
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B. Model simulations 

1. Point charges in random potential. 

The Eqs. Q27I - 130|) provide the practical scheme to 
calculate the optical response of the system in presence of 
the random potential. We begin by considering the case 
of the usual 2D charged Bose gas in random potential. 

Two especially interesting cases of the impurity po- 
tentials have been examined. The first one is the case 
of Coulomb impurities, located at a certain distance h 
over the Cu02 planes: u(r) = Uq (R 2 + ft. 2 ) -1 / 2 , where 
R is the in-plane projection of r. The 2D Fourier trans- 
form of the potential is u(q) = Uo exp(—qh)/q. The 
second case we considered here corresponds to the De- 
bay impurities, u(r) = Uq exp(— r/rO)/r, where rg is 
the screening length. The Fourier transform is u(q) = 

twi + ? 2 r 2 r l/2 . 

The overall view of dimcnsionless spectral dependen- 
cies of optical conductivity, effective relaxation rate T(co) 
and the effective plasma frequency fl p (oj) 003) is shown 
in Fig. [3J It can be seen, that the conductivity spec- 
trum is characterized by a broad asymmetric feature due 
to plasma oscillations of quasiparticles in the potential 
wells at the plasma frequency, and is not very sensitive 
to the choice of the profile of the wells, specified by the 
impurity potential. The stronger is the localization, the 
sharper is the plasma peak. Because of the total spectral 
weight conservation, J Q u(lu) dtv = const, this implies 
the transfer of the intensity from the coherent response 
of the condensate oc 5(to) to the plasma modes. 



To be noted is a strong frequency dependence of effec- 
tive relaxation rate T(co) with nearly constant slope in 
a wide range, reminiscent of the marginal Fermi liquid 
behaviour of the cuprates. 

The comparison with the previously studied 3D Bose 
gas^, 3D and 2D Fermi gas^Lifii makes it possible to 
conclude, that the response of localized carriers always 
remains essentially the same. Still more nontrivial is, 
however, the resemblance of these spectra with those, 
obtained for pinned charge density waves^i and Wigner 
crystals^ 2 ". It seems, that there is a certain dualism 
between the collective modes of localized charges and 
pinned phase degrees of freedom. 



2. Skyrmion-like quasiparticles in a random potential. 

As far as we treat the real-space electronic inhomo- 
geneity in oxides within the effective model of 2D Bose 
gas, the basic results of previous section hold also for this 
case. 

However, the topological texture of EHBL phase is 
not the system of structureless point quasiparticles, but 
rather skyrmion-like entities with a complicated internal 
structure. Hence, the explicit form of the effective ran- 
dom potential does not remain as simple as what follows 
from the Coulomb law. Nevertheless, among the possi- 
ble interactions of the charged impurity with the inho- 
mogeneity the direct electrostatic one is of prime impor- 
tance. 

It is possible to model the charge distribution within 
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the quasiparticle in a direct manner, assigning the charge 
density p, to its point (Fig. ??), if the discrete nature of 
the quasiparticle in the CuC>2 is taken into account. The 
interaction of an impurity with the quasiparticle is ob- 
tained, summing over "internal" coordinates of the quasi- 
particle (see Eq. ffity): 

u(r) - Y^piUoKr-n), (31) 

i 

whence the Fourier transform is straightforwardly ob- 
tained: 

u(k) = U f(k)F(k), (32) 

where 

F(k) =J2pi exp(ik-r0 (33) 
l 

I 



where r.y = |r$ — Tj\ is the distance between two arbi- 
trary internal points of the quasiparticle, Jo is the Bessel 
function of the first kind. 

The form-factor of the quasiparticle enables us to in- 
vestigate a number of the charge configurations of par- 
ticular interest, such as the charged linear stripes etc., of 
course with reserve, that only the details of radial charge 
distribution enter the average form-factor. 



It is therefore instructive to consider first the simple 
case, when the charge distribution in a quasiparticle has 
the form of an uniform ring with inner radius R\ and 
outer radius i?2- We denote their ratio as a — R1/R2. 
Changing a continuously from to 1, we can model both 
the disk, the thin hook and the intermediate ring config- 
urations. The continuous counterpart of Eq. (|34|l reads: 



(\F(k)\ 2 )(q) = (27T J rp(r)J (qr)dr\ , (35) 
^ Ri ' 

where p = Q/7r(i?| — R\) - is the (constant) charge 
density, Q is the full charge of the ring. 



A series of model optical conductivity spectra, cal- 
culated for the ring geometry of charge distribution in 
skyrmion-like quasiparticle, is depicted at Fig. 01 It 
shows the spectra of the bosonic quasiparticles (mass 
m = 2m e , charge 2e) distributed within the CuC>2 plane 



is the form-factor of the droplet and f(k) = e h /k for 
Coulomb impurities, randomly spread at a distance h 
over the CuC-2 planes. 



For the applications, we have to calculate the ex- 
pression (\F(k)\ 2 ), averaged over the orientations of the 
quasiparticles. To this end, the additional physical as- 
sumptions about the arrangement of the quasiparticlet 
within the CuC>2 plane is required. Since there are a 
priori no reasons to expect some ordering of the quasi- 
particle, we consider below the case of they random ori- 
entations. In two dimensions, we obtain: 



(34) 

I 

with the concentration n-iD = 3.5 10~ 3 A -2 and trapped 
by the potential induced by Coulomb impurities ran- 
domly located at the height ~ 3A above CuC>2 plane 
with the mean wells depth ~ 1500 cm -1 . The value of 
the static dielectric constant e — 40 was used in the cal- 
culation. The outer radius of rings was fixed to be 30 

A. 

As the aperture of the ring increases, its outer radius 
kept constant, the intensity of the plasma peak progres- 
sively drops because the charge gradually moves away 
from the impurity. In other words, the larger is the 
quasiparticle, the weaker it is trapped by the potential. 
However, the increase of the quasiparticle extent does 
not simply reduce to its gradual derealization. It can 
be seen, that spectrum of the plasma peak of the ring 
decays into the oscillating structure, that becomes more 
pronounced as the ring gets thinner. This means, that 
the inhomogeneous charge distribution in the quasiparti- 
cle leads to the complicated profile of the collective mode 
spectrum with a fine resonant structure. To get insight 
into the nature of these oscillations, we remember, that 
the frequency dependence of the memory function M(uS) 
is obtained, according to Eq. I|30|l . as a result of the inte- 
gration of a complicated density fluctuation correlation 
functions over the momentum space. Taking account of 
the form-factor of quasiparticle can be interpreted as the 
evaluation of the memory function using the same inte- 
grand, as for the case of structureless point particles, now 
with a highly nontrivial fc-dependent plasmon density of 
states, that is tantamount to the modification of the dis- 
persion relation of the quasiparticles as distinct from the 
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FIG. 4: The optical conductivity of the sky rmion- like quasiparticles with the ring-shaped charge distribution in a random 
potential. Left panel illustrates the role of different inner-to-outer radius ratio denoted by a. Right panel illustrates the role of 
different values of the bare relaxation parameter 70 = 0.0, 0.01, and 0.05 eV given a = 0.9 



case of the free Bose gas. 

The result obtained is to be compared with 
the microwave optical conductivity of lightly doped 
La2_ 3; Sr 3 ;Cu04 43 (x = 0.015). The multi-peak structure 
of calculated spectra falls in nearly the same spectral 
range, as that, observed in experiment. The authors of 
the paper |43| regard the resonances in optical conduc- 
tivity to be the manifestation of the pinned Wigner crys- 
tal response of the carriers in CuC>2 planes. We believe, 
that our results could provide another (quantitative) in- 
terpretation of these experimental findings in terms of the 
plasma oscillations of the quasiparticles phase, trapped 
in the roughness of a complicated non-uniform potential 
landscape. 

The memory function approach considered here en- 
ables us to model the optical response of the quasipar- 
ticles due to their collective plasma modes. Of course, 
there are also another channels of the dissipation in a 
phase-separated systems, related to the excitations of 
other degrees of freedom, that are beyond the scope of 
the present analysis. We may try to get idea of the role 
of this additional absorption mechanisms phenomenolog- 
ically, inserting an additional relaxation term in the de- 
nominator of the generalized Drude formula. This trick 
reflects the fact, that different mechanisms of current re- 
laxation act nearly independently. Having in mind to 
minimize the number of model assumptions, we consider 
the following trial formula 



er(cj) = Re ■ 



M(cj) 



170 



(36) 



where 70 = const. Note that the trivial memory function 
M = i"/o corresponds to the conventional Drude formula. 



The right hand side panel in Fig. 0] illustrates the 
effect of a gradual increase of the bare relaxation param- 
eter 70 on the optical conductivity of the "ring-shaped" 
quasiparticles. In absence of 70 the optical conductivity 
displays the broad plasma feature and pseudogap depres- 
sion of the spectrum at lower frequencies. Turning on 70 
leads to the development of the Drude like tail within the 
gap. 

It is interesting to note that such a spectral behavior is 
actually observed for different cuprates (see e.g. Refill) 
with varying the temperature. It may be seen, that the 
inclusion of additional dissipation channel allows to sim- 
ulate the temperature effect. 



V. CONCLUSION 

We consider the cuprates to be unconventional systems 
which are unstable with regard to a self-trapping of the 
low-energy charge transfer excitons with a nucleation of 
EH droplets being actually the system of coupled elec- 
tron CUO4 - and hole CuO^ - centers having been glued 
in lattice due to strong electron-lattice polarization ef- 
fects. The hole/electron doping into parent cuprate is 
likely to be a driving force for a growth of primary EH 
droplets with a gradual stabilization of a single, or multi- 
center skyrmion-like EH Bose liquid collective mode, fol- 
lowed by a first order phase transition from parent in- 
sulating state into unconventional topological EH Bose 
liquid phase. The latter is believed to be described as 
a system of interacting skyrmion-like entities which con- 
centration is specified by the doping. Nanoscopic elec- 
tron inhomogeneity is believed to be inherent property 
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of doped cuprates throughout the phase diagram begin- 
ning from EH droplets in insulating parent system and 
ending by a topological phase separation in EHBL phase. 

We have examined the effects of electron inhomogene- 
ity accompanying such a phase separation on IR optical 
conductivity. A simple model of metal-insulator com- 
posite and effective medium theory has been used to de- 
scribe the static phase separation effects. Specific effects 
of spectral weight red-shift with the appearance of low- 
energy effective insulating gap, the seeming "quenching" 
of effective plasmon frequency, and the emergence of MIR 
bands have been assigned to the effects of phase separa- 
tion. 

The low-frequency dynamics of topological EHBL 
phase in a random potential in underdoped regime has 
been discussed in a quasiparticle approximation in frames 
of the memory function formalism. The appropriate low- 



frequency optical conductivity appears to have anoma- 
lous many-peak "non-Drude" spectral dependence that 
resembles that of observed experimentally. 

The effects of static and dynamic nanoscopic phase 
separation are believed to describe the main peculiari- 
ties of the optical response of doped cuprates in a wide 
spectral range. 
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